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The  aim  of  this  work  is  to  investigate  which  would  be  a  good  preliminary  plantwide  control  structure 
for  the  process  of  Hydrogen  production  from  bioethanol  to  be  used  in  a  proton  exchange  membrane 
(PEM)  accounting  only  steady-state  information.  The  objective  is  to  keep  the  process  under  optimal  oper¬ 
ation  point,  that  is  doing  energy  integration  to  achieve  the  maximum  efficiency.  Ethanol,  produced  from 
renewable  feedstocks,  feeds  a  fuel  processor  investigated  for  steam  reforming,  followed  by  high-  and  low- 
temperature  shift  reactors  and  preferential  oxidation,  which  are  coupled  to  a  polymeric  fuel  cell.  Applying 
steady-state  simulation  techniques  and  using  thermodynamic  models  the  performance  of  the  complete 
system  with  two  different  control  structures  have  been  evaluated  for  the  most  typical  perturbations.  A 
sensitivity  analysis  for  the  key  process  variables  together  with  the  rigorous  operability  requirements  for 
the  fuel  cell  are  taking  into  account  for  defining  acceptable  plantwide  control  structure.  This  is  the  first 
work  showing  an  alternative  control  structure  applied  to  this  kind  of  process. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  current  industrial  hydrogen  production  technology  could 
only  partially  meet  the  requirement  of  the  small-scale  fuel  pro¬ 
cessors,  which  must  be  compact,  turnkey  and  of  high  efficiency 
by  combining  component  functionalities  and  eliminating  unneces¬ 
sary  components.  Although  the  on-board  fuel  processor  program 
has  been  put  on  hold  by  the  U.S.  Department  of  Energy  (DOE), 
a  number  of  new  approaches  for  the  establishment  of  hydrogen 
refueling  station  based  on  some  form  of  a  fuel  processor  are  being 
pursued  [1].  Particularly  among  all  the  alcohols,  methanol  is  the 
most  popular  fuel  for  reforming,  not  only  because  it  requires  mild 
reforming  conditions  and  has  potential  for  attainment  of  highest 
possible  efficiency,  but  also  because  of  its  possibility  to  be  pro¬ 
duced  from  renewable  resources  [2,3].  Another  alcohol  is  ethanol, 
which  is  already  used  as  one  of  the  main  additives  for  gasoline 
for  both  its  high  octane  number  (RON  and  MON)  and  low  toxicity, 
has  gained  popularity  for  its  environmentally  friendliness  as  well 
as  sustainability  [4].  Biomass,  including  bioethanol  and  other  high 
molecular  weight  materials  such  as  sugar  alcohol  is  considered  to 
be  promising  raw  material  in  terms  of  their  carbon  cycling  neutral¬ 
ity.  Conventionally,  hydrogen  can  be  obtained  by  high  temperature 
gasification/pyrolysis  of  biomass  followed  by  catalytic  reforming 
of  the  gas/liquid  product  [5].  A  good  recent  review  about  different 
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technologies  for  integrated  fuel  processors  for  fuel  cell  application 
is  given  in  Qi  et  al.  [6].  They  emphasized  that  process  intensification 
technologies  such  as  engineered  catalysts,  on-site  heat  produc¬ 
tion/removal  and  product  purification  can  not  only  allow  precise 
control  of  reaction  and  heat/mass  transfer  rates,  but  also  help  opti¬ 
mize  the  operation  conditions,  and,  consequently,  improve  overall 
efficiency  and  mitigate  the  requirement  for  materials  and  capital 
investment. 

On  the  other  side,  feedback  control  of  fuel  cell  power  systems 
has  recently  started  to  attract  attention.  Several  control  problems 
for  fuel  cell-powered  electric  vehicles  are  outlined  in  Powers  and 
Nicastri  [7].  Boettner  et  al.  [8]  have  identified  control  opportuni¬ 
ties  for  the  compressor  within  the  fuel  cell  system.  Control  of  the 
electrical  power  output  and  oxidant  supply  in  electrical  power  gen¬ 
eration  for  a  fuel  cell  powered  vehicle  is  discussed  in  Mufford  and 
Strasky  [9].  Pukrushpan  et  al.  [10]  have  derived  a  lumped  dynamic 
model  of  the  cell  stack,  and  regulated  the  net  power  output  by  con¬ 
trolling  the  air  supply  to  the  cathode.  More  recently,  a  comparison 
of  PI  and  LQG  controls  for  air  supply  is  presented  in  Rodatz  et  al. 
[11].  In  Pukrushpan  et  al.  [12]  have  demonstrated  a  model-based 
multivariable  control  design  for  the  fuel  processing  system  (FPS)  to 
regulate  the  temperature  of  CPO  and  the  mole  fraction  of  hydrogen 
in  the  anode.  Gorgun  [13]  presented  a  control-oriented  physics- 
based  modeling  of  FPS  reactors  in  fuel  cell  power  systems.  In  both 
last  references  methane  is  used  as  the  feed  and  dynamic  models 
were  employed  for  defining  the  proposed  control  structure. 

As  the  authors  understand,  up  to  now,  there  is  no  works  con¬ 
sidering  the  overall  FPS  with  PEM,  using  bioethanol,  where  the 
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Nomenclature 

£rev  reversible  open  circuit  voltage 

/  molar  flow  rate 

F  Faraday  constant 

FA  water  flow  to  the  reformer 

FE  ethanol  flow  to  the  reformer 

FEC  ethanol  flow  to  the  burner 

h  specific  heat 

zCeli  load  current  cell 

^Cell  actual  electrical  power  generated 

^system  net  power  produced  by  the  fuel  cell  system 

Qpem  heat  produced  by  the  cell 

R  water/ethanol  molar  ratio 

Rg  universal  gas  constant 

T# 23  temperature  of  stream  #23 

U  manipulated  variables 

VCeU  actual  cell  voltage 

W  perturbations 

Wcomp  work  of  the  compressor 

XE  ethanol  molar  fraction 

Y  controlled  process  variables 

Z  performance  variables 

Acronyms 

DMC  dynamic  matrix  control 

ESR  ethanol  steam  reforming 

FPS  fuel  processor  system 

HHV  higher  heating  value 

FITS  high  temperature  water  gas  shift 

LFIV  lower  heating  value 

LTS  low  temperature  water  gas  shift 

PEM-FCS  proton  exchange  membrane— fuel  cell  stack 

PrOx  preferential  oxidation 

WGS  water  gas  shift 

Greek  letters 
r\  efficiency 

cp  drop  potential  due  to  the  operation  irreversibility 

Subscripts 
a  anode 

Burn  burner,  post-combustion  reactor 

c  cathode 

Cell  fuel  cell  stack 

Et  ethanol 

Ref  reformer  reactor 

Rev  reversible 

Superscripts 

in  inlet  condition 

out  outlet  condition 


main  control  structure  is  analyzed  through  a  steady-state  model.  It 
is  remarked  because  generally  the  steady-state  model  is  the  only 
available  information  for  doing  preliminary  tests.  Flence,  in  this 
work,  a  control-oriented  steady-state  model  of  an  ethanol  proces¬ 
sor  for  hydrogen  production  and  fuel  cell  operation  is  developed. 
The  rigorous  model  consists  of  mass  and  energy  balances,  chem¬ 
ical  equilibrium  and  feasible  heat  transfer  conditions  elements. 
The  FPS  converts  a  bioethanol  fuel  to  a  hydrogen  (H2 )  rich  mix¬ 
ture  that  is  directly  fed  to  the  proton  exchange  membrane  fuel  cell 
(PEM-FC).  Since  ethanol  can  be  produced  from  renewable  feed¬ 
stocks  or  agriculture  residues,  is  an  attractive  option  as  a  fuel 


source.  The  overall  plant  studied  here  includes  an  ethanol  steam 
reforming  (ESR),  water  gas  shift  (WGS)  and  preferential  oxidation 
(PrOx)  reactors  for  the  fuel  processing,  followed  by  a  PEM-FC  and 
a  post-combustion  unit.  The  conceptual  design  of  this  plant  was 
presented  by  Francesconi  et  al.  [14]  in  order  to  achieve  a  highly 
heat  integrated  system.  The  study  presented  there,  about  cost  and 
performance  requirements  of  the  total  power  plant,  was  extremely 
useful  in  understanding  the  system  level  interactions,  implications 
on  system  performance  and  for  control  structure  synthesis  since 
the  conceptual  design  stage.  Therefore,  a  proper  plantwide  control 
structure  design  is  proposed  to  keep  the  optimal  operation  point. 

In  this  paper,  some  control  theoretic  tools  such  as  sensitivity 
analysis  to  critical  perturbations,  as  given  by  Luyben  and  Luyben 
[15]  and  Luyben  et  al.  [16],  as  well  as  some  assumptions  suggested 
in  Pukrushpan  et  al.  [12]  are  accounted.  Even  though  in  this  last 
work  the  feed  to  the  FPS  is  methane  and  did  not  consider  any  heat 
integration  system  some  analogous  aspects  were  helpful  as  a  guide 
to  understand  better  both  steady-state  and  dynamic  behavior.  In 
this  context,  the  problem  is  how  to  handle  the  integrated  system 
efficiency  accounting  the  best  operating  conditions  of  the  reformer, 
the  effect  of  the  inlet  temperature  to  the  WGS  and  PrOx  reactors, 
combustion  preheating,  expander  unit,  crude  ethanol  as  fuel  and 
the  PEM  specifications. 

The  main  operating  variables  were  determined  by  Francesconi 
et  al.  [14]  considering  that  the  endothermic  nature  of  the  reformer 
has  a  significant  effect  on  the  overall  system  efficiency.  The  high¬ 
est  energy  consumption  is  demanded  by  the  reforming  reactor,  the 
evaporator  and  re-heater  operations.  To  obtain  an  efficient  integra¬ 
tion,  the  heat  exchanged  between  the  reformer  outgoing  streams 
of  higher  thermal  level  (reforming  and  combustion  gases)  and  the 
feed  stream  should  be  maximized.  Another  process  variable  that 
affects  the  process  efficiency  is  the  water-to-fuel  ratio  fed  to  the 
reformer.  Large  amounts  of  water  involve  large  heat  exchangers  and 
the  associated  heat  losses.  A  net  electric  efficiency  around  32%  was 
calculated  based  on  the  ethanol  higher  heating  value  (FIHV).  The 
responsibilities  for  the  remaining  68%  are:  dissipation  as  heat  in 
the  PEM-FC  cooling  system  (38%),  energy  in  the  fuel  gases  (10%) 
and  irreversibility  in  compression  and  expansion  of  gases  (20%). 

In  this  work,  a  commercial  simulator  was  used  to  solve  the 
mass  and  energy  balances,  and  to  compute  the  operating  condi¬ 
tions  for  the  process  units  under  the  most  critical  disturbances.  The 
process  under  study  was  implemented  within  FIYSYS®  flow  sheet¬ 
ing  software  [17].  Maximum  heat  integration  within  the  system 
is  necessary  to  achieve  acceptable  net  electrical  efficiency  levels. 
The  use  of  process  integration  techniques  applying  pinch  analy¬ 
sis  has  already  been  reported  for  the  design  of  a  heat  exchange 
network  for  integrated  PEM  fuel  cell  systems  [18].  In  this  paper, 
the  heat  exchanger  network  involving  the  reformer,  burner,  gas 
cleaning  units  and  the  PEM-FC  stack  was  modeled  using  the  LNG 
exchanger  model  which  is  a  tool  available  in  FIYSYS®.  This  allows 
identifying  the  most  successful  heat  exchange  opportunities,  and 
to  define  the  optimal  operating  conditions  of  the  ethanol  processor 
for  obtaining  the  best  global  efficiency  considering  the  plant  bal¬ 
ances.  The  analysis  considered  here  would  be  useful  as  a  basis  for 
further  process  dynamic  simulation  working  at  the  optimal  struc¬ 
ture  determination  and  testing  whether  or  not  the  proposed  control 
structure  works  well  for  each  time  instant. 

2.  Description  of  the  fuel  processor  system 

The  fuel  processing  plant  has  been  built  according  to  Francesconi 
et  al.  [  14]  recommendations.  The  fuel  processor  chemically  converts 
ethanol  to  hydrogen  including  a  hydrogen  cleanup  equipment,  the 
fuel  cell  stack  electrochemically  converts  the  hydrogen  energy  to 
electric  power.  Additionally,  associated  equipment  for  heat,  air  and 
water  management,  and  auxiliary  equipment  such  as  pumps  and 
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blowers  are  accounted.  Fig.  1  shows  the  main  components:  ethanol 
steam  reformer  (ESR),  high  (HTS)  and  low  (LTS)  temperature  water 
gas  shift  reactors,  a  preferential  oxidation  reactor  of  CO  (PrOx),  a 
PEM-FC,  a  combustor,  compressors  and  an  expander.  For  the  steady- 
state  model  the  pressure  drops  are  neglected  and  the  operating 
pressure  is  fixed  at  3  atm.  The  inlet  flows  depicted  in  Fig.  1  are  the 
following:  Stream  #01  is  the  ethanol  flow  provided  at  25  °C.  Stream 
#02  is  the  water  flow  required  for  the  steam  reformer  whose  flow 


Table  2 

Fuel  processor  streams  connected  to  the  LNG  unit. 


rate  is  controlled  by  the  water  to  ethanol  molar  ratio.  Stream  #30 
is  the  airflow  required  for  the  operation  of  combustor,  fuel  cell  and 
the  PrOx  units.  In  the  integrated  system,  the  air  compressor  will  be 
driven  by  a  turbine  expanding  the  post  combustion  gases.  The  com¬ 
pressor,  expander  and  pumps  isentropic  efficiencies  are  75%.  The 
inlet  air  conditions  are  1  atm,  25  °C.  The  most  important  operating 
conditions  are  summarized  in  Tables  1  and  2. 

The  ideal  unit  cell  voltage  is  calculated  through  ( 1 ),  adopted  from 
Godat  and  Marechal  [19]: 


c  _  pO  (t _  \  i  ^Cell 

£  rev  —  %evUCellJ  ^ - 2f -  ’ 


'n(PH2,a)  +  n  *n(Po2,c) 


(1) 


where  Erev  is  the  voltage  calculated  for  the  cell  temperature  (rCell), 
and  a  and  p ^  c  are  the  partial  pressure  for  H2  and  02  averaged 
(arithmetic  mean)  between  the  inlet  and  outlet  conditions. 

Useful  work  (electrical  energy)  is  obtained  from  a  fuel  cell 
only  when  a  current  is  drawn,  but  the  actual  cell  voltage  (VCeii) 
is  decreased  from  its  equilibrium  thermodynamic  potential  (Erev) 
because  of  irreversible  losses.  When  current  flows,  a  deviation  from 
the  thermodynamic  potential  occurs  corresponding  to  the  electri¬ 
cal  work  performed  by  the  cell.  Therefore,  the  expression  of  the 
voltage  of  a  single  cell  is 


^Cell  —  ^rev  —  (2) 

where  cp  is  the  drop  potential  due  to  the  operation  irreversibility 
which  represents  the  deviation  from  the  equilibrium  value.  The 
actual  electrical  power  generated  by  the  cell  (Pceii)  can  then  be 
calculated  from 


^Cell  =  ^Cell  x  icell-  (3) 

The  operating  voltage  is  defined  as  the  power  level  at  which  unit 
cell  voltage  drops  to  0.5  V  from  the  ideal  voltage  [5].  The  current 
z’ceii  is  related  with  the  hydrogen  molar  flow  rate  at  the  anode. 

'Cell  =  2F.Cftn2,a-/H°2Uta)  (4) 


Table  1 

Heating  and  cooling  demands  for  the  system  streams. 


Stream  number 

#04 

#05 

#07 

#09 

#12 

#14 

#28 

Temperature  (°C) 

709 

709 

538.9 

236.7 

406 

80 

80 

Pressure  (atm) 

3 

3 

3 

3 

3 

3 

3 

Mole  flow  (kmol  h_1 ) 

0.03670 

0.06284 

0.06284 

0.06284 

0.06284 

0.06360 

0.17497 

Molar  fraction 

Hydrogen 

0.8000 

0.48602 

0.52206 

0.59209 

0.54734 

0.56642 

0.04117 

H20 

0.0000 

0.28037 

0.24433 

0.17430 

0.18431 

0.15588 

0.15589 

Methane 

0.0000 

0.02561 

0.02561 

0.02561 

0.02445 

0.02523 

0.00919 

CO 

0.0000 

0.11235 

0.07631 

0.06278 

0.00002 

0.00002 

0.00001 

co2 

0.0000 

0.09566 

0.13169 

0.20173 

0.19856 

0.20547 

0.07467 

Ethanol 

0.2000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

Oxygen 

0.0000 

0.00000 

0.00000 

0.00000 

0.00008 

0.00008 

0.08233 

Nitrogen 

0.0000 

0.00000 

0.00000 

0.00000 

0.04525 

0.04682 

0.63667 
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Table  3 

Optimal  values  for  the  key  system  variables. 


Description 

Optimal  value 

^Ref 

Temperature  reactor  ESR 

709  °C 

Fhts 

Temperature  reactor  HTS 

500°C 

Flts 

Temperature  reactor  LTS 

150°C 

FprOx 

Temperature  reactor  PrOx 

237  °C 

R 

Molar  ratio  H2  O-ethanol 

4 

02/CO 

Feed  relationship  02-C0  for 
the  PrOx 

2 

o2/h2 

Flor  relationship  02-H2  at 
the  PEM 

4.76 

F* Cell 

Power  cell 

0.8027 1<W 

Xco 

Molar  Fraction  de  CO 

10  ppm 

Pesr 

Pressure  at  ESR 

304  kPa 

Phts 

Pressure  at  HTS 

304  kPa 

Plts 

Pressure  at  LTS 

304  kPa 

PprOx 

Pressure  at  PrOx 

304  kPa 

^Burner 

Pressure  at  the  burner 

304  kPa 

Fcell 

Cell  temperature 

80°C 

F#Cathode  exhausted  _  F#Anode  exhausted 

Pressure  difference  at  PEM 

Min 

WMem 

Membrane  humidity  at 

PEM 

80% 

By  doing  an  energy  balance  over  the  cell  from  the  inlet  to  the  outlet 
conditions,  the  heat  produced  by  the  cell  which  has  to  be  removed 
is  given  by 

inlets  outlets 

QpEM  =  'y~]fjhj(Tce  11 )  “  Cell)  “  ^Cell  (5) 

i=l  o=l 

It  is  assumed  isothermal  and  isobar  PEM. 

2  A.  Main  process  variables  to  be  accounted 

From  the  work  of  Francesconi  et  al.  [14]  can  be  extracted  which 
are  the  main  process  variables  and  their  optimal  values  since  the 
minimun  energy  consumption  point  of  view  at  the  conceptual 
design  stage.  This  information  is  listed  here  in  Table  3. 

3.  Control  problem  formulation 

As  previously  discussed,  one  of  the  key  requirements  of  the  FPS 
controller  is  to  quickly  replenish  the  hydrogen  that  is  consumed  in 
the  fuel  cell  anode  during  current  (load)  changes.  On  the  other  hand, 
the  FPS  controller  needs  to  reduce  the  H2  generation  when  there  is 
a  step-down  in  the  current  drawn  from  the  fuel  cell,  so  H2  do  not 
to  be  wasted.  This  H2  on  demand  operation  involves  the  following 
objectives  (i)  to  protect  the  stack  from  damage  due  to  H2  starvation, 
(ii)  to  protect  ESR  from  overheating  and  (iii)  to  keep  overall  system 
efficiency  high,  which  includes  high  stack  H2  utilization  and  high 
FPS  ethanol-to-H2  conversion.  Objectives  (i)  and  (ii)  are  important 
during  transient  operations  while  objective  (iii)  can  be  viewed  as 
a  steady-state  goal.  Objectives  (ii)  and  (iii)  are  also  related  since 
keeping  the  desired  ESR  temperature  during  steady-state  implies 
proper  regulation  of  the  oxygen-to-carbon  ratio  which  corresponds 
to  high  FPS  conversion  efficiency. 

The  stack  current,  /load,  is  considered  as  an  exogenous  input  that 
is  measured.  The  control  problem  is  formulated  using  the  general 
control  configuration  shown  in  Fig.  2. 

The  PEM  control  is  formulated  based  on  Feroldi  et  al.  [20]  con¬ 
siderations  based  on  the  advantages  of  using  a  regulating  valve  for 
the  cathode  outlet  flow  in  combination  with  the  compressor  motor 
voltage  as  manipulated  variables  in  a  fuel  cell  system.  At  a  given 
load  current,  the  cathode  input  and  output  flow  rate  determine  the 
cathode  pressure  and  stoichiometry,  and  consequently  determine 
the  oxygen  partial  pressure,  the  generated  voltage  and  the  com¬ 
pressor  power  consumption.  Therefore,  in  order  to  maintain  a  high 


Fig.  2.  Generalized  block  diagram  for  the  plant  control  structure. 


efficiency  during  operation,  the  cathode  output  regulating  valve  has 
to  be  adjusted  to  the  operating  conditions,  specially  marked  by  the 
current  drawn  from  the  stack.  In  addition,  in  [6]  is  demonstrated 
that  the  appropriate  valve  manipulation  produces  an  improvement 
in  the  system  transient  response.  The  implementation  was  done 
with  a  predictive  control  strategy  based  on  dynamic  matrix  con¬ 
trol  (DMC),  using  the  compressor  voltage  and  the  cathode  output 
regulating  valve  as  manipulated  variables.  Hence,  both  the  fuel  cell 
voltage  and  oxygen  excess  ratio  in  the  cathode,  are  regulated  and 
thus,  the  system  performance  can  be  improved. 

3 A.  Assumptions  for  control  problem  simplification 

Several  assumptions  are  made  in  order  to  simplify  the  control  of 
the  FPS.  Since  the  control  of  the  three  cleaning  reactors  (two  water 
gas  shift,  WGS,  of  high  and  low  temperature  and  the  preferential 
oxidation,  PrOx)  are  not  studied,  the  three  components  are  lumped 
together  as  one  volume.  It  is  also  assumed  that  all  of  them  are  per¬ 
fectly  controlled  such  that  the  desired  values  of  the  reactants  are 
supplied  to  the  reactors.  It  is  assumed  that  the  composition  of  the 
air  entering  to  the  blower  is  constant.  Additionally,  any  temperature 
other  than  that  of  the  ESR  is  assumed  to  be  determined  optimally 
by  the  LNG  HYSYS  tool  which  represents  a  perfectly  controlled  heat 
exchange  network.  Finally,  all  gases  obey  the  ideal  gas  law  and  are 
perfect  mixtures. 

Based  on  all  of  these  considerations,  it  is  assumed  as  the  most 
critical  process  variables  to  be  controlled  are  the  ESR  temperature, 
the  relationship  of  H2 O-ethanol  at  the  feed  plant  and  the  feed  flow 
rate  of  H2  to  the  PEM.  In  Fig.  2  can  be  seen  the  generalized  block 
diagram  representing  the  main  plant  control  variables  accounted: 

W  (perturbations):  feed  molar  fraction  (H2 O-ethanol),  current 
demands  and  H2  feed  flow  rate; 

U  (manipulated  variables):  the  pre-heated  ethanol  flow  rate  valve, 
the  ethanol  feed  flow  rate  to  the  burner  and  the  water  recycle  feed 
flow  rate  from  the  PEM; 

Z  (performance  variables):  TRef,  R  and  feed  flow  rate  to  the  PEM; 

Y:  controlled  process  variables. 

Therefore  the  control  objectives  are  to  keep  the  performance 
variables  under  optimal  values,  given  by  the  synthesis  stage, 
through  variations  on  the  manipulated  variables  even  though  per¬ 
turbations  occur. 

3.2.  Proposed  plant  wide  control  structures 

Based  on  the  previous  study  in  the  following  sections  two  con¬ 
trol  structures  are  analyzed  based  on  a  sensitivity  evaluation  to  the 
most  critical  perturbations. 
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Fig.  3.  Control  Structure  No.  1. 


3.2 A.  Structure  No.  t 

This  structure  can  be  seen  in  Fig.  3  which  consists  on  controlling 
molar  relationship  between  ethanol-water  through  a  ratio  con¬ 
trol,  the  ESR  temperature  through  the  heat  duty  given  from  the 
exhausted  gases  of  the  burner.  The  ethanol  feedflow  rate  to  the  ESR 
is  manipulated  according  to  the  H2  requirements  for  the  PEM. 

3.2.2.  Structure  No.  2 

This  structure  can  be  seen  in  Fig.  4  which  consists  on  controlling 
molar  relationship  between  ethanol-water  through  a  ratio  control, 
the  ESR  temperature  through  ethanol  feed  flow  rate  to  the  burner, 
and  its  temperature  through  the  heat  duty  given  by  the  exhausted 
gases  to  the  ESR  (T# 23).  The  ethanol  feedflow  to  the  ESR  is  manipu¬ 
lated  according  to  the  H2  requirements  for  the  PEM. 

4.  Results 

A  chemical  plant  typically  has  a  large  number  of  units  with  mul¬ 
tiple  recycle  streams.  Many  different  control  strategies  are  possible, 
and  it  would  be  impractical  to  perform  a  detailed  dynamic  study  for 
each  alternative.  An  alternative  procedure  is  to  screen  out  poor  con¬ 
trol  structures  which  can  be  done  accounting  with  a  steady-state 
process  model  trying  to  deal  with  problems  such  as  snowball  effect, 


Fig.  4.  Control  Structure  No.  2. 

produced  sometimes  by  material  recycle.  The  principal  objective  is 
to  use  a  steady-state  analysis  that  can  reveal  structures  that  require 
large  changes  in  manipulated  variables  when  load  disturbances 
occur  or  when  a  change  throughput  is  made.  These  structures  can  be 
eliminated  from  further  study.  The  idea  is  to  specify  a  control  struc¬ 
ture  (fix  the  variables  that  are  held  constant  in  the  control  scheme) 
and  specify  a  disturbance.  Then,  solve  the  nonlinear  algebraic  equa¬ 
tions  to  determine  the  values  of  all  variables  at  the  new  steady-state 
condition.  For  realistically  complex  processes,  analytical  solution 
is  out  of  the  question  and  numerical  methods  must  be  used.  The 
use  of  modern  software  tools  (such  as  SPEEDUP,  EIYSYS®,  or  GAMS) 
which  make  these  calculations  relatively  easy  to  perform  is  rec¬ 
ommended  by  Luyben  [21].  Elence,  in  order  to  justify  the  control 
structures  selection  the  FPS  modeled  in  EIYSYS®  is  tested. 

Follow  this  methodology  and  assuming  as  critical  disturbances 
the  water  molar  fraction  in  the  feed  flow  rate  and  the  current  load 
demand  for  both  structures  are  analyzed.  As  can  be  seen  from 
Tables  4-7  how  these  perturbation  impact  on  the  manipulated  vari¬ 
ables  FA,  FE  and  FEC  and  on  the  efficiency.  Here  the  overall  efficiency 
of  the  FCS  is  the  same  defined  by  Francesconi  et  al.  [14]  as  the  net 
energy  output  of  the  system  obtained  from  the  gross  output  by 
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Table  4 

Sensitivity  analysis  for  control  Structure  No.  1. 


Base  case 

XE  (—1%) 

Variation  (%) 

XE  (—5%) 

Variation  (%) 

XE (—10%) 

Variation  (%) 

Disturbance  in  water  molar  fraction 

XE  (ethanol  molar  fraction  in  %) 

1 

0.99 

-1.000 

0.95 

-5.000 

0.9 

-10.000 

FA  (kmolh-1) 

2.936E-02 

2.955E-02 

0.647 

3.039E-02 

3.508 

3.156E-02 

7.493 

FE  (kmolh-1) 

7.340E-03 

7.388E— 03 

0.654 

7.597E-03 

3.501 

7.890E— 03 

7.493 

FEC  (kmolh-1) 

7.581  E-04 

7.878E-04 

3.918 

9.159E-04 

20.815 

1.086E-03 

43.253 

Efficiency 

0.2881 

0.2881 

0.000 

0.2880 

-0.035 

0.2878 

-0.104 

Table  5 

Sensitivity  analysis  for  control  Structure  No.  2. 


Base  case 

XE  (—1%) 

Variation  (%) 

XE  (—5%) 

Variation  (%) 

XE  (—10%) 

Variation  (%) 

Disturbance  in  water  molar  fraction 

XE  (ethanol  molar  fraction  in  %) 

1 

0.99 

-1.000 

0.95 

-5.000 

0.9 

-10.000 

FA  (kmolh-1) 

2.936E-02 

2.968E-02 

1.090 

3.094E-02 

5.381 

3.220E-02 

9.673 

FE  (kmolh-1) 

7.340E-03 

7.420E-03 

1.090 

7.734E-03 

5.368 

8.049E— 03 

9.659 

FEC  (kmolh-1) 

7.566E— 04 

7.566E— 04 

0.000 

7.566E-04 

0.000 

9.060E— 04 

19.746 

Efficiency 

0.2881 

0.2880 

-0.035 

0.2888 

0.243 

0.2885 

0.139 

Table  6 

Sensitivity  analysis  for  control  Structure  No.  1. 


Base  case 

P(-l%) 

Variation  (%) 

P(-5%) 

Variation  (%) 

P(-10%) 

Variation  (%) 

Disturbance  in  current  load 

XE  (ethanol  molar  fraction  in  %) 

1 

1 

0.000 

1 

0.000 

1 

0.000 

FA  (kmolh-1) 

2.936E-02 

2.906E— 02 

-1.022 

2.788E— 02 

-5.041 

2.642E-02 

-10.014 

FE  (kmolh-1) 

7.340E-03 

7.265E— 03 

-1.022 

6.970E-03 

-5.041 

6.605E— 03 

-10.014 

FEC  (kmolh-1) 

7.581E-04 

7.511  E-04 

-0.923 

7.222E-04 

-4.736 

6.875E-04 

-9.313 

Efficiency 

0.2881 

2.880E-01 

-0.035 

2.880E— 01 

-0.035 

2.879E-01 

-0.069 

Table  7 

Sensitivity  analysis  for  control  Structure  No.  2. 


Base  case 

P(-l%) 

Variation  (%) 

P(-5%) 

Variation  (%) 

P(-10%) 

Variation  (%) 

Disturbance  in  current  load 

XE  (ethanol  molar  fraction  in  %) 

1 

1 

0.000 

1 

0.000 

1 

0.000 

FA  (kmolh-1) 

2.936E-02 

2.905E-02 

-1.056 

2.778E-02 

-5.381 

2.618E-02 

-10.831 

FE  (kmolh-1) 

7.340E-03 

7.263E-03 

-1.049 

6.945E-03 

-5.381 

6.545E-03 

-10.831 

FEC  (kmolh-1) 

7.581E-04 

7.566E-04 

0.000 

7.566E-04 

0.000 

7.566E— 04 

0.000 

Efficiency 

0.2881 

2.879E— 01 

-0.069 

2.877E— 01 

-0.139 

2.875E-01 

-0.208 

subtracting  the  electrical  energy  needed  to  operate  FCS  auxiliaries 
such  as  pumps  and  compressors  divided  by  the  heating  value  of  the 
ethanol  consumed  in  the  fuel  processor  for  reforming  and  burning, 


HHV  ^system 

n??s  ~  HHVEt.  x  (FE  +  FEC) 


(6) 


The  HV  factor  can  be  the  higher  heating  value  (HHV)  or  the  lower 
heating  value  (LHV)  of  the  ethanol.  HHV  represents  the  actual 
amount  of  chemical  energy  in  the  fuel  (relative  to  standard  con¬ 
ditions),  while  LHV  neglects  heat  below  150  °C.  In  Tables  2-5  the 
HHV  are  included. 

From  the  simulation  results  can  be  concluded  that  Structure  No. 
1  does  not  present  significant  variation  on  the  efficiency  however 
the  porcentual  variations  on  the  manipulated  variables  are  more 
important  than  those  produced  using  Structure  No.  2.  In  this  last  one 
the  efficiency  is  kept  closed  to  the  same  value  of  the  base  case.  The 
stronger  effect  is  detected  for  the  FEC  when  a  disturbance  at  molar 
fraction  of  about  10%  occurs,  the  FEC  reaches  the  40%  of  variation.  It 
indicates  that  Structure  No.  1  is  the  most  sensitive  since  actuators 
saturation  are  probable  to  occur  during  the  operation  plant  and, 
as  a  consequence,  a  lost  of  controllability.  For  Structure  No.  2  the 
variation  is  less  than  that  produced  in  Structure  No.  1  even  though 
its  magnitude  indicates  to  be  careful  on  the  controller  to  be  chosen 
for  this  specific  case. 


5.  Conclusions 

An  efficient  ethanol  processor  depends  on  the  ability  to  keep 
the  process  at  operating  conditions  of  the  reformer  and  their  effi¬ 
cient  energetic  integration  obtained  from  the  conceptual  design. 
The  influence  of  some  critical  perturbations  can  produce  that  sev¬ 
eral  variables  move  away  from  the  desired  point.  Hence,  accounting 
with  a  good  steady-state  model  helps  for  screening  several  control 
structures  by  a  disturbances  sensitivity  analysis.  This  methodol¬ 
ogy  is  useful  for  showing  the  capability  of  the  model  to  evaluate 
alternative  control  strategies.  Hence  this  preliminary  exploration, 
accounting  the  main  factors  of  the  most  critical  elements  of  the 
plant,  is  useful  for  reducing  the  control  dimensionality  which  can 
be  further  investigated  if  a  control-oriented  dynamic  model  can 
be  accounted  for  the  later  test.  As  it  is  understood,  this  is  the  first 
work  which  proposes  a  plantwide  control  for  this  kind  of  process. 
In  future  works  it  is  thought  to  include  the  dynamic  behavior  of  the 
other  auxiliary  equipment  and  perform  a  deeper  study  on  the  final 
control  structure  and  the  controller  synthesis. 
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